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Abstract 

The rates of electron scattering via phonons in the armchair single-wall carbon nanotubes 

■ ^ , were calculated by using the improved scattering theory within the tight-binding approximation. 

,^ \ Therefore, the problem connected with the discrepancy of the scattering rates calculated in 

c/3 ' the framework of the classical scattering theory and ones predicted by experimental data was 

■ • clarified. Then these results were used for the solving of the kinetic Boltzmann equation to 

describe electron transport properties of the nanotubes. The equation was solved numerically 

C^ , by using both the finite difference approach and the Monte Carlo simulation procedure. 
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1. Introduction 

At present the interest to the study of carbon nanotubes increases dramatically [1-23]. 
CN \ First of all it is caused by prompt development of nanoelectronics technologies. Due to 

^ \ these achievements it is possible to fabricate the pure carbon nanotubes with high-quality 

O ' structure [18,19]. The prospect of their future application as active and passive elements 

2^ ' of nanoelectronics is out of doubts today [6,8,11-17] . Moreover, it is necessary to note that 

the theoretical study of carbon nanotubes presently has much to be desired. It concerns 
both the descriptions of the physical properties of nanotubes, and the descriptions of the 
charge carrier transport in such structures. In our opinion, this problem is connected on 
"Ti \ the one hand by that the nanotubes are specific objects of quantum nature. That is why 

the study of them is rather complicated task. And on the other hand these structures 
O ■ have been investigated intensively only one decade. 

One of the important areas in investigations of the nanotubes is the study of their elec- 
k> \ trophysical properties, especially the properties of the armchair single-wall nanotubes. It 

^ \ is known that such nanotubes can be reproduced in the best way, and their electrophysical 

properties are practically independent on their diameter [14,15,21]. This cannot be said 
about the nanotubes with another chiralities. So, for example, zigzag type nanotubes can 
have either metallic or semi-conducting properties depending on their diameter [1] . That 
is why we will consider armchair nanotubes only. Therefore, to avoid the account of the 
edge effects let us consider rather long structures {L > 20 nm), i.e. we will consider the 
armchair nanotubes with the length greater then electron mean free path [6, 14,21]. 
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To describe the electron transport in such carbon nanotubes the semiclassical approach 
and the kinetic Boltzmann equation for one-dimensional electron gas can be used [6,24]. 
Both the finite-difference method [6,24] and the Monte Carlo imitative simulation [25-28] 
can be used for solving of this equation. Nevertheless, the correct use of these approaches 
is possible only when all the dominating charge carrier scattering processes are described 
in the most precision way. In single-wall nanotubes these scattering processes are the 
phonon ones [3-6,10,14,20-23]. At the same time, as far as we know, this problem is very 
poorly studied yet. Hitherto, in particular, the reason of a huge difference in theoretical 
calculations and experimental definition of electron-phonon scattering rates as well as 
the independence of the latter from the diameter of nanotubes [14, 21] are not cleared 
up. That is why the purposes of the present paper are to calculate the rate of electron 
scattering via phonons in the armchair single-wall nanotubes and to develop a model of 
electron transport in such nanotubes based on the solution of the Boltzmann transport 
equation. 

2. Theory 

Let us consider the phonon scattering in nanotubes at the electric quantum limit, i.e. 
when the angular momentum of all electrons is equal to zero [10]. To meet this condition 
the nanotubes with small diameter d < 3 nm are considered. Such nanotubes correspond 
to (n, n) armchair nanotubes with the chirality index n < 20 [1] . In this case the scattering 
rates can be calculated by using the perturbation theory [29]. Then, taking into account 
the one-dimensional nature of the electron gas in nanotubes, the following equations are 
valid [20,24,27,30-33] 
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where k is the electron wave vector; q is the phonon wave vector; h is the Plank constant; 
E is the electron energy relative to Fermi level Ep which is taken as zero {Ep = 0); g 
is the matrix element of electron-phonon coupling; ks is the Boltzmann constant; cu^/^ is 
the cyclic phonon frequency; T^ is the phonon gas temperature; z/ is the index signing the 
phonon branch; e/a are the indexes signing phonon emission or absorption. Furthermore, 
it should be noted that in eq. (pQ) the Pauli principle is not taken into account [29]. The 
latter can be taken into consideration while the Boltzmann equation is being solved or 
while the Monte Carlo simulation of the electron ensemble drift checking whether the final 
state is free [25-27]. Moreover, the formula for the angular momentums of electrons and 
phonons have not appeared in the system of eqs. Q and (J4)). This formula is not consid- 
ered by us because in the electric quantum limit the electrons will be scattered by only 
the phonons with zero angular momentum due to the angular momentum conservation. 
Therefore we will be taking into account only such phonons. 

The relation between the electron wave vector and the electron energy is given by the 
known formula [2,37] 




E^^2{k) = tJo{1-2cos[ — ]], (5) 



where Ei^2 is the electron energy in the band 1 and 2, respectively [2,4,15]; a = 0.246 
nm is the lattice constant [1,5]; Jq = 2.7 eV is the overlap integral [7,9]. The matrix 
element of electron-phonon coupling can be calculated by using the classical scattering 
theory within the tight-binding model [2,34,35], but only with some improvements which 
we are dwelling on below. In this case the basic equations for the calculation of the matrix 
element can be taken from [2,35] 

m = — 7 — rT2COS ^ — ^ sm ^ — ^ , (6) 
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where I and t are the indexes signing the longitudinal and transverse phonon modes, 
respectively [2]. 

Having calculated the scattering rates of all possible phonon scattering processes in the 
armchair nanotubes we have concluded that the dominant scattering mechanisms are the 
longitudinal optical (LO) and acoustic (LA) phonon backscattering processes as well as 
the transverse acoustic (TA) phonon (twiston) backscattering. The longitudinal phonon 
scattering is intraband with the transition from one Dirac point to another [2-4, 10, 15], 
whereas the TA scattering causes interband transition in the vicinity of Dirac points. The 
other phonon scattering processes can be neglected because their rates are very little in 
comparison with the TA, LA and LO backscattering rates of active electrons, the electrons, 
which are in a vicinity of the Dirac points K and K' close to the Fermi level. It should 
be noted that in other studies only these dominant scattering mechanisms are taken into 
account [4,6,14,20-23]. Thus, the scattering rates in the armchair single-wall nanotubes 
can be calculated according to eqs. (HJ - (0) and the following formulae describing the 
dispersion relations between the phonon wave vector q and the phonon energy huj for the 
phonons with zero angular momentum 
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There are C^^ = 0.459 eV, C^^ = 0.645 eV, C^^ = 0.029 eV, C^"^ = 0.359 eV, Q 
4.740 eV, C^^ = 17.28 eV, C^^ = 23.33 eV, C}f = 13.77 eV, C^^ = 3.146 eV, C^S 
0.918 eV, C^S = 1.100 eV, C^S = 0.195 eV, C^"^ = 0.020 eV, C^S = 0.196 eV, C^^ 



0.069 eV, C^^ = 0.059 eV, C[f = 0.094 eV, C^^ = 0.083 eV, Cj^ = 0.124 eV, C- 
0.023 eV and Cq^ = 0.101 eV in eqs. (jHl) - (fTUj) . These equalities have been found by 



using the analytical approximation with high precision for the phonon-zone structure of 
the armchair single-wall nanotubes. This zone structure is determined by means of the 
zone folding method [37] conformably to phonon-zone structure of graphene [36]. The 
validity of such an approach is discussed in [23] . 

Now let us consider the values fj,i, fj,t, qi and qt. Due to our investigations it is found 
out that there are values fii = 2mc/a, and fit = 4mc/a {ttlc = 2 • 10~^^ kg - the mass of 
carbon atom) in l^f^l and \gt\ , respectively, instead of the value of the armchair single- 
wall nanotube linear mass density fi = Anmc/a [2] (the first improvement of the classical 
scattering theory). It is caused by equality of the carbon atom oscillation phases on 
the whole circumference of the cross-section perpendicular to the nanotube axis for the 
phonon modes with zero angular momentum, i.e. by the synchronism of oscillations of 
all the carbon atoms in a half of the primitive cell. In that case it is possible to consider 
all the atoms in the primitive cell like two, for the longitudinal modes, or four, for the 
transverse ones, massive particles bounded with each other 2n or n times, respectively, 
as strong as two carbon atoms. The discussed phenomenon clarifies the fact that the 
electron-phonon scattering rates in the armchair single-wall nanotubes do not depend on 
their diameters [14,21] as the values fii and /if do not, too, in contrast to the value /i. 

As for the values qi and qt-, which are the deformation constants in the directions 
along axis and circumference of the nanotube, respectively, they can be calculated by 
using the plane-wave approach. In contrast to the well-known Slater local-orbital one 
giving the wrong results (see [4]), where the values qi = qt = qo are determined by means 
of analysis of alteration of the overlap integral Jov due to the distance alteration between 
two carbon atoms in graphene [34,35], in the plane- wave approach the values qi and qt can 
be defined by the same way, but at some different assumptions. The difference is that the 
Slater local-orbital model assumes that the electrons are individual particles belonging to 
individual atoms, whereas the plane-wave model assumes that the electrons are de Broglie 
waves filling the whole space in the crystal lattice of graphene. Moreover, it is supposed 
that the phonon displaces only the ion but it does not displace the electron orbit in crystal 
lattice within the Slater local-orbital approach [34,35], whereas within the plane- wave one 
it is supposed that the phonon displaces the whole carbon atom. In an approach like that, 
in fact, it is taken into consideration the "orbital relaxation" (see [4]) when the relaxation 
time is assumed to be equal to zero. 

Then, using the results of [2,34,35,37], the following formulae can be written for 
the unfolded graphite sheet obtained from the armchair single-wall nanotube within the 
plane-wave approach at the tight-binding approximation 
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Eov{R,k) = -JoJl + Acos^{Rkj2)+Acos{RK/2)cos(^V3Rky/2y (16) 

Here r is the two-dimensional radius-vector; k is the two-dimensional wave vector; x 
is the coordinate along the nanotube axis Ox; y is the coordinate along the nanotube 
circumference Oy. In eq. (fT^ the coefficient equal to 2 appeared because of there are 
two carbon atoms in the primitive cell of graphene. According to eqs. (fT^ - p6|) the 
value of qi is equal to 8.643 nm~^, and the value of qt is equal to 9.980 nm^^. Let us 
note that the values of qi and qt calculated in the framework of the plane- wave approach 
(the second improvement of the classical scattering theory) are less then the value of 
go = 25 nm~^ calculated within the Slater local-orbital approach [35]. Our results are 
in a good agreement with the assumption of [4] in which the values of g; = 8 nm~^ and 
qt = 10 nm~^ were considered so as to secure the accordance of the theoretical results 
with the experimental data. 

For example, in fig.l the functions of the electron-phonon scattering rates versus the 
electron wave vector for the armchair nanotubes calculated by using eqs. (^ - (lllj) at 
temperature T^ = 290 K are presented. It can be noted that the LA and LO phonon 
scattering rates are of the same order whereas the TA phonon scattering rate is about ten 
times less than those at T,y = 290 K. 

Let us write down the equations describing the state of the electron gas in the armchair 
nanotube. In general, these equations can be presented as [6, 20, 24-28] 
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Here t is the time; /i,2 = fi,2{k,x,t) is the electron distribution function in the band 1 
and 2, respectively; I is the operator of the electron-phonon interaction with account of 
the Pauli principle; F = F{x,t) is the electric field strength in the nanotube along Ox; 
1*1,2 = Ui^2{k) is the electron group velocity in the band 1 and 2, respectively; n^ = n^i^x, t) 
is the electron linear concentration; Udr = Udr{x, t) is the electron drift velocity; J = J(x, t) 
is the electric current in the nanotube; e is the value of elementary charge. According 
to [27] the operator / effects on the electron distribution function as follows 

// = E (W^/(l - /') - W^7'(l - /)). (23) 




Here the sum is over all possible final states; /' is the electron distribution function in a 
final state; W is the rate of scattering causing the return of electrons to the initial state. 
At the tight-binding approximation the eq. (fTT^ can be written as (see the eq. ^) 



«i,2 = TVsin — . (24) 




So, to study the kinetic processes in the armchair single-wall carbon nanotubes it is 
necessary to solve the kinetic Boltzmann equation (fTTj) using the equations of motion ()18|] 
and (J2l|) . Then the kinetic parameters ul, Udr and J can be calculated according to eqs. 

dUl) -(1221). 

3. Results and discussion 

Let us consider electron transport in the infinite length armchair nanotubes in uniform 
constant electric field. Such a transport is practically equivalent to electron transport in 
rather a long nanotube placed on the ideal metal contacts [6,14,21]. Meanwhile, the 
nanotube can be considered as long if its length is greater or at least the same order as 
the TA phonon limited mean free path of active electrons. In this case the eqs. (II 7j) and 
(1201) ^^6 reduced to 

riL = -. (26) 

a 

The solving of eq. ()25j] allows the kinetic parameters to be calculated. Specifically it may 

be the electric current in the nanotube. 

In fig. 2 and 3 some results of calculations of the electric current in the armchair 
nanotubes at Tta = T = 290 K are presented. Here T is the temperature of the nanotube 
surroundings. These results are obtained by using both the numerical solution of eq. (j^ 
by the finite difference approximation and the imitative Monte Carlo simulation of the 
electron ensemble drift [38-40]. The results coincide and are not resolved in the figure 
scale. 

Moreover, the experimental data are presented in fig. 2 and 3 along with the theoretical 
results obtained by A. Javey et al. [14], who applied the Monte Carlo simulation procedure 
with use of fitting parameters. These figures evidently demonstrate a very good agreement 
between the theory and experiment. Especially it concerns those results where the non- 
equilibrium LA and LO phonons (Tla,lo > T = 290 K) [23] are taken into account. Some 
discrepancy of the results in fig. 3 at low fields is explained by the TA phonon pinning effect 
at some conditions when the TA phonon mode is partly suppressed due to an influence 
of the substrate on oscillations of carbon atoms in the nanotube [6]. 

We did not consider the pinning effect, which is probably difficult to describe theoret- 
ically, but we considered a most important effect like heating of the LA and LO phonon 
gas, i.e. the non-equilibrium LA and LO phonons [23]. As a first order approximation 
this effect can be described by the following obvious formula found from the energy con- 
versation law 

N (o^K, Tla.lo) - N {uk, T) = r]FJ (F, Tla,lo) , (27) 

where ujk is the cyclic frequency of the LA and LO phonons at K point {Huk = 0.151 eV) 
[36] ; 7] is the parameter which is determined by the quality of the thermal contact between 
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the nanotube and its surroundings. The value of this parameter obviously depends on the 
type of matter of the surroundings and the area of the contact as well as another factors. 

In our opinion, it is a very difficult task to calculate the value of rj theoretically, but it 
is easy to find it from the experimental data by using the fitting procedure to theoretical 
results. So, we did like this. And the results presented in fig.2 and 3 have been obtained 
at such values of t] which give rise to the best agreement between the theoretical results 
and experimental data. In fig.4, in addition, the functions of the LA and LO phonon 
gas temperature versus the electric field strength are given. These functions have been 
calculated by means of eq. (|7fj) at the best fitted parameter rj. It is clear from this figure 
that the non-equilibrium phonons cause the superheating of the LA and LO phonon gas 
(see [23]). 

We would like to add that the TA phonon pinning effect can be described, in principle, 
by the same way like the superheating effect of the LA and LO phonon gas in nanotubes. 
In other words it is possible to consider the pinning effect as if the supercooling of the TA 
phonon gas takes place (Tta < T). Then to achieve the agreement between theory and 
experiment in fig. 3, like it is in fig.2, it must be guessed that the temperature of the TA 
phonon gas Tta is equal to 200 K (T = 290 K). 

So, the validity of the improved scattering theory is verified by the things discussed 
above and the results presented in fig.2 and 3. Moreover, it is confirmed by the following 
reasons. The TA phonon limited mean free path of active electrons It as well as the LA 
and LO phonon limited one /; calculated theoretically are in a good accordance with the 
same values, respectively, found from the experiments. For instance, the experiments, 
which the pining effect is not observed in, predict that the value of k is equal to 500 ± 200 
nm at T = 290 K [14,22]. The improved scattering theory predicts that the value of 
If = Up iW^ji^{kp) + Wj^j^i^kp))' [21] is equal to 632 nm. Here up and kp are the absolute 
value of the group velocity and electron wave vector, respectively, at the Fermi level 
in the armchair single-wall nanotube. It also follows from the experiments that li is 
equal to 10-15 nm at very high electric fields (~ 10'^ V/m) [6,14,21]. The theory gives 
that the quantity of k = up iWlj^{kp) + Wljs^{kp) + WlQ{kp) + Wloikp))'^ possesses 
the following values: Ii{Tj^a,lo = 290 K) = 23.8 nm, //(Tla,lo = 1200 K) = 15.0 nm and 
h{ThA,LO = 1970 K) = 10.0 nm. According to the results discussed above and taking into 
consideration the superheating of the LA and LO phonon gas, one can draw a conclusion 
that the improved scattering theory appropriately describes the real scattering processes 
in the armchair single-wall carbon nanotubes. Besides, it predicts that the temperature 
of the LA and LO phonon gas is between 1000 K and 2000 K at very high electric fields. 
These temperatures are much less and more reliable than those which were guessed in [23] 
(> 6000 K). 

Thus, in the present paper the rates of the electron scattering via phonons in the 
armchair single-wall carbon nanotubes are calculated by using the improved scattering 
theory. The model of charge carrier transport in such nanotubes based on the solving 
of the Boltzmann transport equation is developed. It is found out that the discrepancy 
of theoretical calculations in the framework of the classical scattering theory and exper- 
imental estimations of the values of k and k has the following reasons: (i) the features 
of carbon atom oscillations for the one-dimensional phonon modes with zero angular mo- 
mentum are not taken into account; (ii) the excessive values of qi and qt, obtained by 
Pietronero et al. for graphene [35] by using the Slater local-orbital approximation [34], 
are used; (iii) the superheating of the phonon gas is not taken into consideration. Ev- 
erything noted above sets conditions for that there is a discrepancy between the value of 



li (It) calculated in the framework of the classical scattering theory and the value of // 
(It) obtained experimentally. As an example, this discrepancy is more than two times for 
(10, 10) nanotube. Therefore, to provide the agreement of theory with experiment it is 
necessary to use the values ni and Ht instead of the value fi (the first improvement of the 
classical scattering theory) as well as the value of qi = 8.643 nm~^ and qt = 9.980 nm~^ 
instead of the value of go = 25 nm~^ (the second improvement of the classical scattering 
theory) while calculating the phonon scattering rates of electrons in armchair single-wall 
carbon nanotubes. Moreover, it is also necessary to take into account the non-equilibrium 
LA and LO phonons at high electric fields. 

We think that we would succeed much more in the agreement of the theoretical results 
with the experimental ones if we considered the finite length nanotubes instead of the 
infinite ones. But in that case the computational complexity would increase much more, 
too. So, in conclusion we would like to say that our subsequent purpose is the calculation 
of the electrical characteristics of the armchair single-wall carbon nanotubes at various 
electric fields F = F{t) in the framework of the developed scattering theory. 

The authors gratefully wish to acknowledge Dr. D. Weiss for his careful examination 
of the paper. 
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Figure captions 



Fig. 1. LA, LO and TA phonon scattering rates in the band 1 (a) and 2 (b) of the armchair 
single-wall carbon nanotubes. 

Fig. 2. The function of the electric current versus the electric field strength in armchair single- 
wall carbon nanotubes {tj = 1.063 • 10~^ W~^m). 
The solid curve is the experimental data from [14] (L = 700 nm). 
The circles are the theoretical results from [14] (L = 700 nm). 
The points are our theoretical results at Tla,lo = 290 K (L = oo). 

The squares are our theoretical results in which it is taken into account the non-equilibrium LA 
and LO phonons {L = oo). 

Fig. 3. The function of the electric current versus the electric field strength in armchair single- 
wall carbon nanotubes (?? = 6.737 • 10~^ W~^m). 
The solid curve is the experimental data from [21] (L = 1000 nm). 
The points are our theoretical results at Tla,lo = 290 K (L = cx)). 

The squares are our theoretical results in which it is taken into account the non-equilibrium LA 
and LO phonons [L = oo). 

Fig. 4. The function of the temperature of the LA and LO phonon gas versus the electric field 
strength in the infinite length armchair single-wall carbon nanotubes at 77 = 1.063 • 10~^ W~^m 
(solid curve) and rj = 6.737 • 10^^ W^^m (dashed curve). 
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